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ABSTRACT 


Previous  calculations  of  barium  cloud  prowth  failed  to 
pive  reasonable  estimates  of  striation  scale  size  and  onset 
time  because  they  either  assumed  no  ionospheric  end-shortinp 
or  no  ion  diffusion  (all  modes  unstable)  or  stronp  shortinp 
(all  modes  stable).  We  derive  a  set  of  equations  which  in¬ 
clude  finite  end-shortinp  and  solve  hy  a  new  numerical  tech¬ 
nique.  For  moderate  end-shortinp  (in  the  experimental  ranee) 
we  obtain  a  fastest  prowinp  mode  of  reasonable  size  on  the 
rear  of  the  cloud.  The  prowth  rates  and  instability  ranpe 
are  very  sensitive  to  the  end-shortinp  and  mapnitude  of  the 
ambient  electric  field.  There  is  no  prowth  below  a  critical 
value  of  E  or  above  a  critical  value  of  the  shortine.  A 
series  of  praphs  are  presented.  We  also  derive  most  previous 
results  as  special  cases  of  our  equations. 


m 


TABLF  OF  CONTENTS 


Abstract . 

*  •  •  •  * 

Fist  of  Illustrations  ... 

•  •  •  •  • 

Introduction . 

Basic  Equations  . 

*••••• 

Dynamical  Equations  for  Equilibrium  and 
I erturbations  .  . 

Numerical  Solutions  and  Discussion. 

Results*  A:  St,ecial  Cases  Yiol«<inc  Previous 

**•••• 

Appendix  B:  Computational  Procedure. 

Acknowledgment . 

References.  .  .  .  , 


7 

11 


17 

21 

22 

22 


Preceding  page  blank ' 


v 


LIST  OF  ILLUSTRATIONS 


FiR.  1 
Fig.  2 
Fig.  3 
Fig .  4 
Fig.  5 

Fig.  6 
Fig.  7 

Fig .  8 

Fig.  9 
Fig.  10 
Fig.  11 
Fig.  12 


Equilibrium  growth  for  shorting  ratios  X=0.25, 
4  and  50.  * 


Time  evolution  of  r.m.s.  density  perturbation  for 
various  k*  (Qt=10,  X»4,  H*=l,  F*=-10). 

Growth  of  k*=4  perturbation  in  x-space  and  time 

Cfii-10,  X  =  4,  F*  =  l ,  F  *= - 1 0)  . 

x  y 

Growth  rate  y  (at  t*=0.02)  versus  k*  for  various 
shorting  ratios  X  (Qt=10,  F*  =  l,  F*=-10). 


Growth  rate  y  (at  t*=0.2)  versus  k*  ^or  various 
magnitudes  of  the  ambient  electric  field  (£Xt=10 
X=4,  F*/F*=10). 

Maximum  growth  rate  y  (at  t*=0.2)  as  a  function 
ot  F*  for  various  valUeS  of  y  ((Qt)=10,  F*/F*=-10'l  . 

Growth  rate  y  (at  t*=0.2)  versus  k*  for  various 
angles  between  F.  and  the  x-direction  (Qt=10,  X=4 
E.*(.'01m/21. 


'A 


Growth  rate  y  (at  t*-0.2)  versus  k*  for  various 
angles  between  FA  and  the  x-direction  and  for 
|F-aI  =  5.0  and  3.0  respecti velv  fnT'in,  X  =  4)  . 


Growth  rate  y  (at  t*»0.02)  versus  k*  for  various 
values  of  Hr  C A  =  4 ,  F*  =  l,  F*  =  -10'). 

^  y 

Maximum  growth  rate  yMAX  times  Q/ Tl  +  TiCXt')  2  as  a 
function  of  (Qt)  (X=4,  R*=l,  F*=-101. 

Perturbed  density  contours  for  the  k*=4  mode  at 
t*=0.2  CJlT*10,  X-4,  F*  =  1 ,  F*  =  -101. 

Perturbed  potential  (gj)  contours  for  the  k*=4 
mode  at  t*=0.2  (Qt=10,  X  =  4,  EJ  =  1,  f*  =  - 101 


Preceding  page  blank 


vi  i 


There  seems  to  be  general  apreement  that  the  pxR  insta- 

ww*  *w 

hility  ^Simon,  196^]  (sometimes  called  the  pradient  drift 
instability)  is  responsible  for  the  onset  of  striations  in 
barium  clouds  released  in  the  ionosphere.  However,  the  problem 
has  been  to  demonstrate  that  this  is  actually  so  by  calculatinp 
the  scale  size  and  delay  time  and  makinp  comparison  with 
experiments.  Previous  attempts  failed  because  they  either 
predict  instability  for  all  wavelenpths  and  with  comparable 
prowth  rate  [Linsow  and  Workman, 1970 ;  Volk  and  Haerendel .  1971] 
or  because  they  predict  stable  behavior  fSimon.  19701. 

We  have  previously  derived  a  set  of  equations  fSimon  and 
Sleeper,  1972]  which  povern  the  prowth  of  a  barium  cloud 
immersed  in  a  conductinp  medium.  We  now  show  that  all  earlier 
results  are  included  as  special  cases  of  these  equations  cor- 
respondinp  to  no  end  shortinp  or  stronp  shortinp  or  neplect 
of  diffusion.  In  each  of  these  cases,  one  can  account  for  the 
results  obtained.  We  have  now  perfected  a  new  numerical  scheme 
for  solvinp  these  equations  with  arbitrary  end-shortinp .  We 
find  that  there  is  a  preferred  wave  number  and  prowth  rate  and 
that  the  results  are  very  sensitive  to  the  amount  of  end 
shortinp  and  to  the  strenpth  o^  the  external  ambient  electric 
field.  For  values  comparable  to  those  in  actual  releases,  we 
find  a  prowth  rate  and  scale  size  in  peneral  apreement  with 
the  observations. 

Tn  the  next  section,  we  rederive  the  peneral  equations 
from  a  sliphtly  more  peneral  viewpoint.  Followinp  that  we 
describe  a  particular  solution  of  these  equations.  The 
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corresponding  numerical  solutions  and  a  discussion  of  the 
results  are  then  given.  Finally,  in  Appendix  A,  we  show  that 
previous  results  are  indeed  special  cases  of  our  set  of  equa¬ 
tions,  while  Appendix  B  gives  details  of  the  computational 
procedure . 


BASIC  EQUATIONS 

We  shall  denote  the  number  density  of  barium  ions,  ambient 
ions  and  electrons,  respectively,  by  nR,  nA  and  nQ.  Each 
obeys  the  continuity  equation 

+  v-  (n  v)  =  o  m 

and  quasineutrality  holds,  i.e. 

I  nB  +  r,A_nel  «•  ne  (2) 

where  we  assume  singly  charged  ions,  for  simplicity.  Denoting 
electrical  current  density  por  species  as  .TR  ■  e  nRvR ,  JA  = 
e  nAv A,  Je  *  -e  neVe  and  using  Eqs.  Cl)  and  C2),  one  has 

V'  ifA  +  Jb  +  =  0  (3) 

and  this  equation  can  be  used  in  place  of  the  electron  con¬ 
tinuity  equation. 

The  velocities  are  given  by  the  usual  diffusion  relations 
in  a  magnetized  medium 

nv„  =  ±/tnE„  co 

C-^n±/4njgi3  +  csj 
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written  in  the  rest  frame  of  the  neutral  background.  The 
upper  signs  refer  to  ions  and  the  lower  to  electrons  and  all 
symbols  have  their  usual  meaning  in  the  neutral  rest  frame. 
The  corresponding  species  current  divergence  is 


±  %[-D7„n  ±^ne„  ) 

*  £  J.J  i  JJX t,ar)  vn  *  tj.  b 


(*) 


where  we  neglect  temperature  variation  in  the  nerpendicul ar 
direction  (over  the  cloud  dimensions)  and  assume  an  electro¬ 
static  field. 

If  one  now  substitutes  F.q.  (6)  in  Eq.  (3),  and  divides  all 
terms  by  De,  one  obtains  a  complicated  expression  which,  however, 
takes  on  a  simple  form  in  the  limit  (At)  -►».  Tn  this  limit, 

Eq.  (3)  becomes 

V„  Un£  +  -frnc£„  J  =  o  (7, 

'e 

We  assume  that  moving  along  a  field  line  passing  through  the 
cloud,  one  ultimately  reaches  a  non- conducting  region  (the 
lower  atmosphere)  where  all  currents  must  vanish.  Hence,  we 
integrate  Eq.  (7)  twice  and  obtain 

ne  -  <i>  =  -  $<rL)  an 

where  $  is  the  electrostatic  potential  and  g  is  a  function  of 
rj_only.  This  result  for  #  is  of  course  only  the  first  term 
in  a  systematic  expansion  of  the  equations  in  powers  of  T)  /n 
[Simon  and  Sleeper,  1972], 


-4- 


Now  integrate  Fq .  (3)  along  a  magnetic  field  passing 
through  the  cloud  all  the  way  from  one  non-conducting  terminus 
to  the  other.  By  the  boundary  condition  this  becomes 

J  Vj-'IJia  -+  Zx&  +  Zxe.  3  dJ  =  0 

where  the  subscript  1  denotes  the  last  two  terms  in  Fq.  (6) 
and  dl  is  a  differential  length  along  the  magnetic  field 
direction.  Let  us  further  divide  this  integral  into  that  part 
which  lies  in  the  cloud  and  that  which  lies  in  the  ionosphere 
outside  the  cloud.  Then  we  write 

IdJ -I'JsJj.  [£a  -bj^e  +X,e  ]  d%  —  0  ^ 

C  load  ionO 

To  lowest  order  in  I)+/r>e,  the  species  expression  for  VI-JI, 
obtained  by  substitution  of  F.q.  (8)  in  Fq.  (6),  is 

-T  -  ±  Vj.-(-VjV±n  3-  +  ^ne) 

j.W)<Lnx  Vjj  b  3  JJj.ir\T)JkvJn*vlr\i  b  nn) 

C  Mg 


We  now  assume  that,  in  the  cloud  portion  of  the  integral  in 
Hq.  (9),  the  ambient  ions  are  negligible  compared  to  the 
barium  ions.  Thus, 


rig  —  Hg 

_  (cloud  region) 

Y)a  S  0 

and  that  the  reverse  is  true  in  the  ionsophere  portion, 


ne 

(ionosphere  region) 

* 

With  this  assumption,  Fq.  (9),  evaluated  to  lowest  order  in 
n+/ne>  takes  the  form: 
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+  ^  t  riu  - -f^  7j^]-+ -£  Aj.h  -  c> 


Mere  we  have  defined 


A  s  -L  nB  d* 


SOx  S  vj  +  =  D?  ( )  +  %■ ) 


X  =•  }Jj 


Cd  =  -  /  (  »;  +  -2  Of  ;  ^  n 

£  =  e  /  (p?+Mf)n,  di 

°  J;Ono 

h#  s  -  §•/„„„  ^ 

and  neglected  variation  of  diffusion  coefficients  over  the 
dimensions  of  the  cloud.  We  have  allowed  for  election  cross1 
field  transport  contributions  in  the  lower  ionosphere  in  the 
definitions  of  ^T,  E  and  A,  Mote  that  E^  is  the  height  inte- 
prated  (actually  along  a  field  line)  Pedersen  conductivitv 
of  the  ionosphere. 

One  obtains  a  complete  set  of  equations  by  combining 
Fq.  (11)  with  the  barium  ion  continuity  equation  integrated 
over  the  length  of  the  cloud  and  evaluated  to  lowest  order 
in  n+/ne, 

^  ^  C  1  f  _  (T\  i  I  li  l  I  S\  "1 


-/Aj.inr)c  r/sf*b  =  o 
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and  the  continuity  equation  for  the  ambient  ions.  These  repre¬ 
sent  a  set  of  coupled  integro-dif ferential  equations  in  the 
variables  Nc(rx,t),  g(i^,t)  and  n^(r,t).  We  now  partially 
decouple  these  by  the  assumption  that  variations  in  the  ambient 
ion  density  have  small  influence  on  the  dynamics  of  the  ion 
cloud.  This,  in  effect,  neglects  the  back  effect  of  image 
striations  [Volk  and  Haerendel,  1971]  on  the  cloud.  With 
this  assumption,  we  set  n^  in  Eq.  (11)  equal  to  its  value  in 
the  undisturbed  ionosphere  which  means  that  V.n,*n,  rir*0, 

V^Ep-0  (all,  over  the  dimensions  of  the  cloud)  and  ER 
itself  is  the  height  integrated  Pedersen  conductivity  of  the 
undisturbed  ionosphere.  Eq .  (11)  reduces  to: 

%  ■  ['A  A/ -fAi  N  f  ]  -f  v7>  k  '  Vx  A/ 


where  we  drop  the  subscript  on  N  and  (Pi)  and  where  this 
equation  plus  Eq.  (13)  comprises  a  complete  set.  The  boundary 
conditions  appropriate  to  an  ion  cloud  created  suddenly  in  a 
localized  region  at  t»0  are 


N  o 

as 

rx 

i 

e* 

1 

*  M* 

as 

^  oo 

(15) 


where  E^  is  the  ambient  electric  field  (in  the  neutral  rest 
frame).  Eq.  (14)  is  essentially  equivalent  to  Eq .  (8)  of 
Simon  and  Sleeper  T19721  while  Eq .  (13)  is  the  integrated 
form  of  Eq.  (1)  of  that  naner.  In  that  paner,  these  equations 
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werc  solved  in  the  limit  Ep>>euj_N,  i.e.  hiphlv  conductinp 
backpround  ionosphere.  In  this  paper,  we  describe  a  new 
numerical  technique  which  allows  us  to  solve  the  equations 
for  aribtrary  values  of 


DYNAMICAL  EQUATIONS  FOR  EQUILIBRIUM  AND  PERTURBATION 
We  shall  solve  Eqs.  (13)  and  (14)  in  two  stapes  (which  are 
entwined  in  the  actual  numerical  procedure) .  First  a  numerical 
solution  for  a  time  and  space  dependent  cloud  eouilibrium  with 
slab  peometry  is  obtained.  If  we  choose  the  B-  field  in  the 
z-direction  and  assume  that  the  density  and  electric  fields 
vary  only  in  the  x-direction,  F.q.  (14)  becomes: 

jl  _ u  NX\  J.  ik  iXlk  _  _4e.  tX  =  o 

3xL^j3X  J  (flT)  e  3X“  " 


Intepratinp  this  once  from  the  interior  of  the  cloud  to  a 
distant  point,  usinp  the  boundary  conditions  in  Eo.  (15), 
one  has 


2l‘-  X-  + 

j+x  int)  J+* 


where 

>v(X,-0 


<Z  U  j.  Cx,*) 

^  8 


(17) 


Note  that  -3p°/3y  ••=■  EAy  by  the  vanishinp  of  VxF .  This  solution 
can  then  be  substituted  in  Eq.  (13)  which  becomes: 


(18) 


0 


where 
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CX(X,t)  —  - - - 

J  -4* \(X,4.) 


(19) 


It  is  a  straightforward  matter  to  solve  Eq.  (18)  numeri¬ 
cally.  Note  that  a  represents  the  effect  of  end-shorting  and 
that  this  varies  from  point  to  point. 

Next,  we  linearize  Eqs.  (13)  and  (14)  about  the  equilib¬ 
rium  derived  above  and  fourier  analyze  in  the  v-direction  (exn  iky) . 
These  equations  are 

- &[  -  4 ‘n,  ]  ■*  ">  41’  3 

+  =  0 


and 


■ "  --  "alii 

1 


Mi-  ^4  W 

(OX  )  3  * 

-HfAii 


1 1 


JjL 

e 

A  I  J  i 


(21) 


We  solve  Fqs.  (20)  and  (21)  by  studying  the  growth  in  time 
of  some  arbitrary  perturbations  in  x.  Details  of  the  numerical 
procedure  are  given  in  Appendix  B. 
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Throughout  our  numerical  procedure  we  assume  that  the 
initial  density  distribution  is  a  gaussian  with  characteristic 
length  h.  Thus 


-'(x/h) 


n  «■*  o)  =  si  e 


(22) 


Using  this  characteristic  length  h  and  the  initial  central 
density  IT  we  define  the  dimensionless  variables 

*'=  ,  i  *=  -t. Qj./h1 

x*  ■=  e^a/S-B  ,  n*=.h/g 

i*=  Si  A  ,  fjuj./jou 

and  5*= 

cF . 

and  transform  to  the  frame  moving  with  the  velocity  — .  The 
dimensionless  form  of  Fqs.  (18)  -  (21)  are 


2N 

dt 


12. 

inr) 1  dx 


A/  )  ~  0 


(27,) 


c< 


_J _ 

/ 


(24) 


(28) 
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4fJi x(H+& -t-irc^x]  ?i  =  o 

where  all  quantities  are  in  dimensionless  units  and  we  have 
suppressed  the  asterisk.  Also 


(26) 


n / 

d* 


£a>  .  /\  ..  itL  i  £a±  Kd. _ 

i-*MI  a/  <>*  >-^aa7 


(27) 


in  the  same  units.  In  practice,  we  work  with  Fq.  (26)  and 
the  difference  of  Eqs.  (25)  and  (26),  which  is 

;  hi) A-  (C*rt?2tL  9  :=  0 

Note  that  the  entire  behavior  of  a  Riven  k  perturbation 
depends  only  on  Xr  FAx,  EA  ,  and  fix  for  the  barium  ions.  Tt 
is  also  worth  notinp  that  the  equilibrium  elective  field  in 
the  cloud  interior  (in  dimensionless  units)  follows  from  Fq. 
(8)  and  is 


(28) 


£x  - 

=  £ 


"Te/Tg 


Ai 


)  -b(r</r6) 


jl  2tL 

a/  T>T 


Note  also  that  the  boundary  conditions  on  n1  and  Rj  are  that 
n^O  and  Vp^O  as  |x|-*». 
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NUMF.RTCAL  SOLUTIONS  AND  DISCUSSION 


Solutions  arc  obtained  by  integrating  forward  in  time  using 
an  initial  gaussian  ^or  the  equilibrium  and  an  initial  density 
perturbation  consisting  of  a  sum  of  several  Uermite  polynomials 
with  amplitudes  appropriate  to  white-noise  and  with  randomlv 
selected  phases.  A  linear  instability  is  detected  by  ob¬ 
serving  the  growth  in  time  of  the  average  density  disturbance 
defined  as: 


'"Ik" 


n  1  k  (x  *  t) 


|nlk(x,o)  |2dx 


1/2 


(29) 


Details  of  the  comrmtational  procedure  arc  given  in  Appendix 

B.  Let  us  now  examine  a  number  of  the  results. 

Variation  with  X .  To  study  the  effect  o^  variations  in 

X,  we  fixed  the  remaining  variables  at  (fir)*10,  Fx"l,  Ev*-10. 

The  magnitudes  of  (0t)  and  |F.|  are  roughlv  in  the  ranee  of 

experimental  releases  while  the  E  ,  E  ratio  is  arbitrary. 

a  y 

Some  typical  time  developments  of  the  equilibrium  solution  are 
shown  in  Fig.  1  for  three  values  of  X;  0.25,  4  and  50  corres¬ 
ponding  to  .  tronp  shorting,  a  typical  experimental  value  and 
weak  shortinp,  respectively.  Note  that  in  the  strong  shorting 
case,  there  is  diffusion  at  the  ion  porpendi cular  rate  as 
expected  [Simon,  1955)  and  that  very  little  distortion  occurs 
since  the  shorting  ratio  a  remains  close  to  unity  at  all  times. 
The  drift  to  the  right  is  due  to  the  combined  effect  of 

2 

Pedersen  ion  current  and  the  actual  nAy  drift  being  (fir)  cE^/ 
Bfl+(flT)2]  rather  than  cO^/D.  In  the  more  typical  second 
case,  we  see  reduced  diffusion  and  distortion  of  the  cloud 
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due  to  the  varying  value  of  a  across  the  cloud.  Finally,  in 
the  weak  shorting  case,  we  see  a  clo.td  controlled  almost 
entirely  by  electrons  (ambipolar  cross-field  diffusion). 

Since  we  assumed  (fiT)e-*-«»,  there  is  little  change  in  time. 

The  time  behavior  of  some  perturbations  about  the  X«4 
equilibrium  are  shown  in  Fig.  2.  Note  that  the  growing  modes 
rapidly  approach  an  exponential  growth  and  that  there  is  a 
fastest  growing  mode  in  the  vicinity  of  k«s4.  Very  short  and 
very  long  wavelengths  are  stable.  In  fact,  in  a  much  longer 
timescale,  the  growth  rates  of  all  the  crowin"  modes  decrease 
and  some  become  stable  as  the  equilibrium  diffuses  and  its 
density  gradient  decreases.  It  is  interesting  to  note  that 
the  growing  modes  are  localized  on  the  rear  of  the  cloud  as 
may  be  seen  from  the  spatial  evolution  of  the  k*4  mode  as 
shown  in  Fig.  3. 

Some  idea  of  the  effect  of  variation  in  X  upon  these 
results  is  seen  from  the  plot  in  Fig.  4  of  the  exponential 
growth  rate  y^, 


"I 

Ya  =  <n<>> 


<J<  "//$> 

- — 

d*t 


130) 


versus  k  from  various  values  of  X  (y^  evaluated  at  t«0.02). 

Note  the  strong  dependence  on  X.  For  no  end-shorting  X*®, 

all  values  of  k  are  unstable  and  y  approaches  an  asymptotic 
2 

value  as  k  -**».  This  is  in  qualitative  agreement  with  the 
results  in  Linson  and  Workman  [1970]  and  indeed  with  the 
result  in  Simon  [1963]  if  one  takes  the  limit  (fix)  in  that 
paper.  As  X  decreases,  the  growth  rage  decreases  and  a  finite 
range  of  unstable  k  is  established.  For  strong  shorting  all 
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modes  are  stable  for  X  slightly  less  than  0.5. 

It  is  interesting  to  compare  the  magnitude  of  the  growth 
rate  with  the  simple  estimate  [Linson  and  Workman,  19701 

y  ^  (3D 

In  dimensionless  units,  this  is 

y 1+  (AT)*  p*  ^  )OD  (32) 

(at) 

in  our  case.  We  see,  from  Pig.  4,  that  the  growth  is  about 
seven  times  this  value  for  and  rapidly  decreases  as  X 
decreases . 


Variation  with  |F.a|.  Next  we  fix  the  angle  between  P-A 
and  the  slab  normal  at  FAy/EAx«- 10 ,  set  X-4  and  vary  I P-A I  •  The 
resulting  growth  rates  (at  t=0.2)  are  shown  in  Pig.  5.  Mote 
the  rapid  decrease  in  both  maximum  growth  and  range  of  unstable 
k  as  | Ea |  decreases.  In  fact,  if  one  plots  versus  Ey,  as 

shown  in  Fig.  6,  one  obtains  a  straight  line  which  crosses  the 
axis  at  about  F*«-0.5  For  X*y.  Note  that  as  X  decreases,  the 
critical  value  of  |PA|  increases.  The  linear  dependence  indi¬ 
cates  that  Y»4Ax  does  indeed  vary  directly  as  c  FA/Rh.  Indeed 
in  this  range  of  X,  varies  roughly  as 


aA  ^ 


(33) 


where  a,  h,  c  and  d  are  constants. 

Variation  with  Electric  Field  Angle.  If  we  define  0= 
tan"’1(Ev/Fx)  ,  the  variation  of  y  (at  t=0.2)  with  0  is  shown 
in  Fig.  7.  The  growth  rate  is  largest  at  angles  close  to  tt/2 
and  -tt/2  fthc  -tt/2  curve  is  identical  to  that  for 
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+v/2].  It  is  somewhat  higher  for  the  positive  values.  411 
modes  are  stable  at  0=0.  This  is  not  a  general  result.  Insta¬ 
bility  occurs  at  6=0  Tor  large  values  of  |F.a|  for  small  (nr). 

A  plot  of  the  angular  variation  for  two  other  values  of  |F.| 
is  shown  in  Fig.  8. 


Variation  with  (fit) .  The  variation  of  y  with  fir  is  shown 
in  Fig.  9.  One  should  remember  however  that  the  dimensionless 
growth  rate  includes  (Qx)  in  its  definition.  In  real  time 

Y  _  Y  *  c(T'.+  Te)  inr; 

e  B~hL  (nr)z 

r-  »  (34) 

_  Y*  c  (nr)_ 

5  max  e*  Bh  i-ttnt)1 

2 

Hence,  we  plot  y^ax(0t) / [1+ (Or)  ]  versus  (Dx)  and  observe 
that  the  real  growth  rate  increases  with  increasing  (Ox)  and 
approaches  a  constant  value  for  large  (fix).  This  is  shown  in 
Fig.  10.  Note  that  this  is  close  to  the  value  c  F^/Rh  (for 
X»4) . 


Shape  of  Fastest  Mode.  We  note  that  the  maximum  value  of 
k  for  values  of  X  in  the  neibhborhood  oT  x=4  are  o^  the  order 
of  k=4.  This  corresponds  to  a  true  wavelength  in  the  v 
direction  which  is  Xy  =  2rh/k*  s  J  h.  Thus,  the  fastest 
growinr  mode  will  have  comparable  dimensions  in  both  the  x 
and  y-  directions.  A  plot  of  the  density  perturbation  is 
given  in  Fig.  11.  The  corresponding  perturbation  in  g^ 
is  shown  in  Fig.  12.  Note  that  the  maxima  are  displaced 
a  quarter-wavelength  from  the  density  maxima. 

Discussion .  In  summary,  we  have  demonstrated  that  one 
must  include  finite  ionospheric  end-shorting  of  the  cloud 
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in  order  to  obtain  a  most  favored  growth  made  of  reasonable 
size  and  that  the  results  are  quite  sensitive  to  the  shorting 
ratio.  Previous  results  can  be  explained  in  terms  of  either 
no  end  shorting  or  no  diffusion  (all  k  unstable")  or  large  end¬ 
shorting  (all  k  stable).  Shorter  wavelengths  are  stabilized  by 
ion  diffusion  at  moderate  values  of  X.  'Ve  demonstrate  that 
the  growth  rate  is  also  a  sensitive  function  of  electrid 
field  strength  and  angle  and  that  the  plasma  is  stable  if 
|E |  drops  below  a  critical  value.  In  addition,  the  growth 
rate  increases  with  increasing  Ht  but  becomes  a  constant 
for  large  values  of  same.  The  fastest  growing  mode,  for 
moderate  values  of  X  has  a  k  about  4  corresponding  to  fairly 
circular  striations  in  the  plane  perpendicular  to  P,  and  grows 
with  a  rate  of  the  order  of  c | | /B  h. 

One  could  relate  this  to  delay  time  by  estimating  the 
number  of  e-folding  times  required  for  a  thermal  fluctuation 
to  grow  and  appreciably  modulate  the  original  cloud.  Py 
elementary  statistical  mechanics,  the  r.m.s.  fluctuation  of 
the  number  of  ions  contained  in  a  striation  of  length  I,  and 
transverse  dimensions  h  is 

A/s/  L  kl  L  Iri1  J 

where  n  is  the  average  ion  number  density  in  the  cloud.  Taking, 
as  crude  estimates,  n  5  If)6  cm"'\  h  *  1  km,  I,  *  10  km,  one 
finds  AN/N  ■  10  ^  •  e  If  we  thus  assume  that.  25  e-folding 

times  are  required  for  a  fluctuation  of  the  fastest  growing 
k-mode  to  reach  observable  size,  we  obtain  an  estimate  of 
the  delay  time  (in  dimensionless  units)  as, 
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(75) 


Using  the  value  of  Y*^y  *  60  for  |F^|=10  and  X-4,  one  would 

estimate  a  delay  time  of  1/2  in  dimensionless  units  or  about 

half  the  time  to  diffuse  across  the  original  cloud  size  since 

1  2 

t  in  real  time  is  t^  =  yh  /J)^ .  This  is  of  the  order  of 

tens  of  minutes  in  many  experiments. 

In  conclusion,  we  note  that  a  two-dimensional  calculation 
would  differ  from  the  one -dimens i onal  one  in  one  important  way. 
The  finite  ends  of  the  cloud  in  the  v-direction  would  then 
allow  polarization  charges  to  accumulate  and  one  would  have 
a  mechanism  for  shielding  F;^v  in  the  interior  with  consequent 
steepening  of  the  cloud  rear.  This  would  add  an  additional 
steepening  mechanism  which  should  also  be  strongly  \ -dependent . 
One  might  expect  the  time  scale  for  this  steepening  to  be  at 
least  as  fast  as  the  cross-field  diffusion,  so  that  striation  and 
steepening  would  he  stronglv  entwined.  Tn  fact,  if  the  steep¬ 
ening  was  rapid  compared  to  diffusion  spreading,  one  would  have 
a  mechanism  which  enhanced  the  instability  growth  rate  as 
time  went  on  and  a  time  plot  such  as  that  given  in  Fig.  2 
might  well  show  the  curves  bending  upwards  instead  of  down. 
Further  insight  awaits  the  results  of  2-D  calculations  which 


are  now  under  studv. 


-17- 


Appendix  A 

Special  Cases  Yielding  Previous  Results 


No  Fnd  Shorting.  We  show  here  that  one  recovers  the 
results  of  Linson  and  Workman  [1070]  in  the  limit  where 
This  limit  corresponds  to  X-*®  and  thus  F,q .  [23)  reduces  to 
3N/3t=0  with  the  general  solution 


M  —  N  (1) 


[A .  1 ) 


The  correspondinr  value  of  c°  mav  he  obtained  from  Pn .  f27). 
In  this  limit 


^  UJ  -  ~L  d N  4- 

t)  x  N  '3x  O.T 


[A.  21 


We  make  corresnondence  with  Pinson  and  Workman  hv  the  choice 


n  c  x )  .=.  e 


[A.  31 


Then 


23^  -  /  -  =  t  .on*  t  Ant 

?  X  (IT 


[A. 41 


and  Fq  .  [26)  becomes: 


ithi  _  r  aL  lA  ( i  _ 

T**  ^  nxKI  at1  J 


[A.. SI 


while  Po .  [28)  becomes 


jn,  IWjWl//  -  id-Vi  l  :  A  t/^1, --M?  =  0 

-r>  af  ^  nr  /yli~r  L  K  (ir  vi  r a . 6 ) 
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Since  these  are  linear  equations  in  n^  and  Ne^ ,  with  constant 
coefficients,  we  fourier  analyze  each  of  these  variahles  as 
exp[iut  +  imx]  and  immediately  ohtain  a  dispersion  relation. 
The  result  is 


id 


A  £. Ay  illlrtT)']  i i  J 

(OT)  l  A  +  M(0T)±L(nL+j2l)(nT) J 


(A. 7) 


which  one  easily  recognizes  as  precisely  the  dispersion  re¬ 
lation  in  Linson  and  Workman  [3970  ,  Pq.  OS')]  correctinp  for 
the  fact  that  the  ahovc  result  is  in  the  frame  movinp  with  the 
cnAy/B  velocity  and  is  in  dimensionless  units.  Thus,  the 
results  of  Pinson  and  Workman  are  for  the  case  of  no  shortinp 
by  the  ambient  ionosphere  which  implies  electron  (amhipolar) 
diffusion  across  the  magnetic  field  (in  their  actual  calcu¬ 
lation  ami  strong  instability  since  the  stabilizing 

effect  of  end-shortinc  is  removed.  It  should  he  noted  that 
finite  electron  diffusion  f(r2t)c  finite  rather  than  infinitel 
does  stabilize  verv  short  wavelengths,  but  this  is  for  lengths 
very  much  shorter  than  those  stabilized  by  ion  diffusion 
(by  the  ratio  f  (Jit")  +/ (Hr) pl  and  much  smaller  than  striation 

lenpths  seen  in  actual  experiments. 

Infinite  Fnd  Shortinp.  In  this  case  and  A+D. 

Prom  Eq .  (26),  we  see  that  in  this  limit 


-  =  o 


(A. 8) 


The  only  solution  that  vanishes  at  x=±”  is 


f.  *  o 


(A.  01 


-19- 


By  the  equilibrium  condition  in  Eq.  (27) ,  we  have 


Lohi-tan-t 


(A. in) 


Hence  Eq.  (25)  becomes: 


-+  (A1-  iAtnr =  o  fA.m 

which  obviously  has  decavinp  solutions  only.  Hence,  all  per¬ 
turbations  are  stable  in  the  case  of  infinite  end-shortinp  (and 
in  the  limit  r)+/r>e*n,  as  was  assumed  in  the  derivation  of  the 
basic  equations). 

Ion  Diffusion  Neglected.  We  recover  the  results  of  Volk 
and  Haerendel  (1971)  i f  we  use  the  approximations  in  their 
paper.  These  are  (in  our  notation), 

A/  n  A/d  exp  K  x ) 

E.AX  =  o  (A. 12) 

0  '  L  ifx  “  7$  ®  0 

and  ncplect  of  n_£,Vn  compared  to  p^nE.  This  last  assumption, 
in  dimensionless  notation,  means 


»  l 

We  now  apply  these  to  Eo.  (25)  and  obtain 


(A.  IS) 


1  ^  ^4y  ^  ^ 


0 


(A. 14) 
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where  we  take  (!>t)»1  and  use  Pq.  (27)  to  see  that  3e”/9x  has 
two  terms  of  order  unity  and  PAy/n-  ,  resnectively.  Now  usinr 
this  result  and  our  approximations  in  Eo.  f 28)  we  obtain 


Mi  -f  f  +  £a*i 

^  L  3X  '  mi+M)  int)  J”,- 


0  (A. IS) 


The  prowth  rate  is  determined  hy  the  Inst  term  onlv.  Usinr 

the  definition  of  the  equilibrium,  we  find  the  rrowth  rate 
Y  to  he: 


Y  =. 


( 


ar)  [/-h-L-i 
L  J 


(A. 16) 


This  is  in  dimensionless  units.  Tn  real  units,  this  becomes 


r  _  -  %  k  ( /  -b  j 


-l 


(A . 17) 


Comparing  this  with  Eq.  (49)  of  Volk  and  Haerendel,  we  see 

that  a  similar  result  obtains  with  Zp/e  y playinr  the  role 

of  their  UR.  This  identification  is  exact  in  the  limit  of 
small  kv. 

Thus,  we  see  that  the  results  in  Volk  and  Uaerendel  do 
include  the  effects  of  end  shotting  However,  their  ncplect 
of  ion  diffusion  removes  this  stabilizinr  effect  on  short 
waveler pths  and  allows  all  these  modes  to  prow  at  the  same 


rate . 
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Appendix  B: 

Computational  Procedure 

Finite  difference  methods  are  used  in  obtaininp  the 
evolution  of  the  equilibrium  density  N  and  the  perturbations  n^ 
and  p^.  Equation  (23)  is  a  nonlinear  parabolic  equation  for 
the  unknown  N.  It  can  be  solved  by  the  Crank-Nicolson  finite 
difference  scheme  with  forward  projection  of  the  nonlinear 
coefficient  [von  Rosenberg,  1969],  The  Crank-Micol son  analogs 
of  (26)  and  (28)  also  provide  two  tridiaponal  systems  for  the 
unknowns  n^  and  The  grid  points  for  evaluating  n^  and  pj 

are  located  at  time  levels  tanAt.  The  prid  noints  for  evalu- 
atinp  N  are  located  at  time  levels  t^in+^OAt,  so  that  the  value 
of  N  together  with  3g°/3x  from  (27)  can  be  easily  used  in 
evaluating  the  coefficient  of  the  two  tri diagonal  systems. 

In  carrying  out  the  computations,  we  first  solve  for 
the  self-consistent  g^  from  the  initial  value  n^.  We  then 
solve  for  N  at  t-^-At  from  its  initial  value.  The  two  tri¬ 
diagonal  systems  arc  then  solved  to  obtain  values  o^  n^  and 
p^  at  t=At.  We  then  update  the  value  of  N  to  t^At  and  solve 
for  n^  and  gj  at  t=2At  and  repeat  these  procedures.  The 
accuracy  of  the  results  was  checked  by  comparing  results 
obtained  using  different  grid  spacings  in  space  and  in  time 
and  using  different  total  number  of  spatial  grid  points  with 
the  same  spacing. 
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